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We have identified several important features in the neutron scattering data of cuprates, which are 
difficult to be explained in terms of d-wave and isotropic s-wave order parameters. Alternatively, 
we show that the neutron data are in quantitative agreement with an order parameter that has an 
extended s-wave (Ai g ) symmetry and opposite sign in the bonding and antibonding electron bands 
formed within the CU2O4 bilayers. The extended s-wave has eight line nodes and change signs when 
a node is crossed. This A\ g pairing symmetry may be compatible with a charge fluctuation mediated 
pairing mechanism. 



The microscopic pairing mechanism responsible for 
high-temperature superconductivity in copper-based per- 
ovskite oxides is still a subject of intense debate 
despite tremendous experimental and theoretical ef- 
forts for over 15 years. The debate has centered 
around the role of antiferromagnetic spin fluctuations 
in high-temperature superconductivity and the sym- 
metry of superconducting condensate (order parame- 
ter). Extensive inelastic neutron scattering experi- 
ments have accumulated a great deal of important data 
that should be sufficient to address these central is- 
sues. Of particular interest is the magnetic resonance 
peak that has been observed in double-layer cuprate su- 
perconductors such as YBa2Cu30 y (YBCO) [1-6] and 
Bi 2 Sr2CaCu208+ a (BSCCO) [7,8], and in a single-layer 
compound Tl 2 Ba 2 Cu0 6+a (Tl-2201) [9]. A number of 
theoretical models [10-13] have been proposed to explain 
the magnetic resonance peak in terms of d-wave magnetic 
pairing mechanisms. These theories can qualitatively ex- 
plain some features of neutron data but are particularly 
difficult to account for an important feature: The mag- 
netic resonance in optimally and overdoped double-layer 
cuprates is much more pronounced in the odd channel 
than in the even channel [2-4]. In order to overcome 
this difficulty, Mazin and co-worker [14] proposed an or- 
der parameter that has isotropic s-wave symmetry and 
opposite sign in the bonding and antibonding electron 
bands formed within the CU2O4 bilayers [14]. However, 
this model predicts that the resonance energy is larger 
than twice the magnitude of the superconducting gap 
along the Cu-0 bonding direction, in disagreement with 
experiment. Further, the nodeless s-wave gap symmetry 
is inconsistent with the measurements of the penetration 
depth, thermal conductivity and specific heat, which con- 
sistently suggest the existence of line nodes in the gap 
function of hole-doped cuprates [15,16]. 

Here we identify several important features in the neu- 
tron scattering data of cuprates. These features are in- 
consistent with d-wave order parameter (OP). Alterna- 
tively, we show that the neutron data are in quantitative 
agreement with an order parameter that has an extended 
s-wave {Aig) symmetry [17,18] and opposite sign in the 
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FIG. 1. a) The imaginary part of dynamic spin suscepti- 
bility as a function of excitation energy for YBa2Cu3C>6.g2 
in the normal and superconducting states. The figure was 
reproduced from Ref. [5]. b) gj-scan at E= 40 meV in 
YBa2Cu3 06.92 displaying a modulation typical of odd ex- 
citation. The figure was reproduced from Ref. [4]. The 
background (lower line and open squares) was obtained from 
q-scans across the magnetic line. The upper solid curve is a 
fit to a + bF 2 {Q) sin 2 (jvzqi), where F(Q) is the Cu magnetic 
form factor. The modulation is not complete as even excita- 
tions are sizable at qi= 3.5, 7 but with a magnitude 5 times 
smaller [4]. 

bonding and antibonding electron bands formed within 
the Cu 2 4 bilayers [14]. 
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The magnetic resonance peak observed in opti- 
mally doped and overdoped double-layer compounds 
YBa 2 Cu 3 O y [1-5] and Bi 2 Sr 2 CaCu20 8+y [7,8] is a sharp 
collective mode that occurs at an energy of 38-43 meV 
and at the two-dimensional wavevector Qaf — (tt /a, 
7r/a), where a is the nearest- neighbor Cu-Cu distance. 
Fig. la shows the imaginary part of the odd chan- 
nel spin susceptibility as a function of excitation en- 
ergy for a slightly underdoped double-layer compound 
YBa 2 Cu306.92- The figure is reproduced from Ref. [5]. 
There are several striking features in the data: (a) The 
resonance peak at E r ~ 41 meV appears below T c and 
this resonance peak intensity in the superconducting 
state lodd(-^r) is larger than the normal-state intensity 
ladder) by a factor of about 3.6, i.e., lf dd {E r )/I? dd {E r ) 
= 3.6; (b) A spin gap feature is seen below T c and there is 
a small shoulder that occurs at an energy slightly above 
the spin gap energy E g =32 meV. Fig. lb shows a gz-scan 
at E= 40 meV in YBa2Cu3C>6.92 displaying a modulation 
typical of odd excitation. The figure is reproduced from 
Ref. [4]. From Fig. lb, we find feature (c): The odd- 
channel magnetic resonance intensity within the CU2O4 
bilayers is larger than the even-channel one by a factor of 
about 5 (Ref. [4]), i.e., I S odd {E r ) / l! ven {E r ) = 5. From the 
neutron studies on different double-layer compounds with 
different doping levels, one identifies feature (d): The res- 
onance energy E r does not increase with doping in the 
overdoped range but is proportional to T c as E r /ksT c ~ 
5.2 in both underdoped and overdoped ranges [8]. 

On the other hand, the neutron data for a single-layer 
Tl-2201 are different from those for double-layer com- 
pounds. For the single-layer Tl-2201, we show in Fig. 2a 
the difference spectrum of the neutron intensities at T = 
27 K (<T C ) and T = 99 K (>T C ), at a wavevector of Q = 
(0.5, 0.5, 12.25). The figure is reproduced from Ref. [9]. 
Although a sharp peak feature is clearly seen at E r = 
46 meV, it is remarkable that the peak intensity in the 
superconducting state is close to the magnetic neutron 
intensity in the normal state, that is, I s (E r ) / I N (E r )~ 
1. This is in sharp contrast with the above result for the 
double-layer compound YBa2Cu30g.92 where we found 
Io dd {E r ) / Ioddi E r) = 3-6. Thus we identify feature (e) as 
I (E r ) I I N (E r ) ~ 1 for the single-layer compound. 

We would like to mention that feature (e) identified for 
the single-layer Tl-2201 is valid only if the nonmagnetic 
background has a negligible temperature dependence be- 
low 100 K. It is known that nuclear contributions pre- 
dominantly constitute the main part of the neutron sig- 
nal [3]. Any nuclear contributions (nonmagnetic back- 
grounds) are only expected to obey the following stan- 
dard temperature dependence, 1/[1— exp(— E/ksT)] (see 
detailed discussions in Ref. [3]). In the energy range of in- 
terest (E > 40 meV), this temperature dependence factor 
is close to unity and nearly independent of temperature 
for T < 150 K. Indeed, the normal-state neutron intensi- 
ties of both YBCO [3] and BSCCO [7,8] show a negligible 
temperature dependence for T c < T < 150 K. For E — 
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FIG. 2. a) The difference spectrum of the neutron inten- 
sities of single-layer Tl-2201 crystals (with a total volume of 
0.11 cm 3 ) at T =27 K (<T C ) and T = 99 K (>T C ), and at 
wavevector Q = (0.5, 0.5, 12.25). The figure is reproduced 
from Ref. [9]. The solid line is guide to the eye. The differ- 
ence spectrum tends to zero at about 56 meV (about 10 meV 
higher than the resonance energy). This suggests that the 
nonmagnetic background at 56 meV has negligible tempera- 
ture dependence below 100 K. b) The difference spectrum of 
the neutron intensities of optimally doped BSCCO crystals 
(with a total volume of 0.06 cm 3 ) at T =10 K (<T C ) and T = 
100 K (>T C ), and at wavevector Q = (0.5, 0.5, -14). The fig- 
ure is reproduced from Ref. [7] . The difference spectrum goes 
to zero at 60 meV (about 16 meV higher than the resonance 
energy). This suggests that the nonmagnetic background at 
60 meV has negligible temperature dependence below 100 K. 

E r = 46 meV in Tl-2201, the factor l/[l-exp(-E/k B T)] 
decreases by 0.4% when the temperature is lowered from 
99 K to 27 K. If we take the nonmagnetic background 
of 2500 counts per two hours [9] , the nonmagnetic back- 
ground decreases by 10 counts per two hours when the 
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temperature goes from 99 K to 27 K, and by 60 counts 
per two hours when the temperature goes from 150 K to 
99 K. The variation of the nonmagnetic background at E 
= 46 meV below 100 K is negligibly small compared with 
the magnetic resonance intensity (^150 counts per two 
hours). Therefore, feature (e) identified for the single- 
layer Tl-2201 is well justified. 

In order to further justify feature (e), we compare the 
ratios of the magnetic to the nonmagnetic background 
signals in those neutron experiments on different com- 
pounds. From the neutron data, we find that the signal- 
to-background ratio is about 6% for Tl-2201 (Ref. [9]), 
about 12% for an optimally doped BSCCO [7], about 
6% for an overdoped BSCCO [8], and about 100% for 
an overdoped YBCO [3]. The signal-background-ratio in 
the overdoped BSCCO is similar to that for Tl-2201. But 
the overdoped BSCCO shows a difference spectrum sim- 
ilar to that for the overdoped YBCO where the signal- 
background-ratio is larger than that for the overdoped 
BSCCO by a factor of 20. This suggests that, if fea- 
ture (e) identified for Tl-2201 were an artifact caused by 
a small signal-background-ratio, one would not have ob- 
served the intrinsic difference spectra for the overdoped 
BSCCO. The fact that the difference spectra for the over- 
doped BSCCO are similar to the one for YBCO sug- 
gests that one indeed finds the intrinsic magnetic differ- 
ence spectra for the overdoped BSCCO even though the 
signal-background-ratio is only 6%. There are no reasons 
to believe that only the difference spectrum for Tl-2201 
is an artifact. Therefore, the pronounced difference be- 
tween the difference spectrum of the single-layer Tl-2201 
and that of the double-layer BSCCO (compare Fig. 2a 
and Fig. 2b) is due to the fact that Io dd { E r) I Io dd (E r ) 
» 1 for the double-layer BSCCO while I s (E r )/I N (E r ) 
<~ 1 for the single-layer Tl-2201. Moreover, we will show 
below that, for the intraband scattering (even channel) 
in YBa 2 Cu 3 6 . 92 , l! ven {E r ) / lg en {E r ) = 0.72, in good 
agreement with feature (e) for the single-layer Tl-2201. 
This consistency gives further support to the thesis that 
feature (e) is intrinsic. 

One may also argue that the glue that glues about 
300 small crystals of Tl-2201 on Al-plates would cause 
a substantial decrease of the nonmagnetic background 
below 100 K. If this argument were relevant, one 
would have observed a similar effect in overdoped 
Y .9Cao.iBa2Cu 3 07_j / (YBCO-Ca) because 60 larger 
crystals of YBCO-Ca are also glued on Al-plates [19]. 
Because the total volume of YBCO-Ca crystals is larger 
than that of Tl-2201 crystals by a factor of 3.2, one can 
readily show that the amount of glue for YBCO-Ca crys- 
tals is comparable with or even larger than that for Tl- 
2201. This suggests that the effect of glue on the nonmag- 
netic background in Tl-2201 is similar to that in YBCO- 
Ca. From the data for YBCO-Ca (Fig. 2a and Fig. 3a of 
Ref. [19]), one can clearly see that the difference spectrum 
at 50 meV is very close to zero. Because the magnetic 
signal at an energy that is about 10 meV higher than the 



resonance energy is independent of temperature below 
100 K (see Fig. la), the nearly zero value of the difference 
spectrum at 50 meV implies that the nonmagnetic back- 
ground at 50 meV has negligible temperature dependence 
below 100 K in the case of YBCO-Ca. Even in the case 
of Tl-2201 (see Fig. 2a ), the difference spectrum tends to 
zero at about 56 meV (about 10 meV higher than the res- 
onance energy) although one needs more data points in 
the vicinity of 56 meV to definitively address this issue. 
This indicates that the nonmagnetic background at 56 
meV has negligible temperature dependence below 100 
K in the case of Tl-2201. 

We would like to mention that the normal-state mag- 
netic intensities in Fig. 4 of Ref. [9] are significantly un- 
derestimated. This is because the authors assume that 
the g-width of the magnetic peak in the normal state is 
the same as that in the superconducting state [9]. This 
assumption is unphysical. From Fig. 19d and Fig. 19h of 
Ref. [6] , one can clearly see that the q- width of the mag- 
netic peak in the normal state is a factor of 1.6 larger 
than that in the superconducting state in the case of un- 
derdoped YBa 2 Cu306.8- Further, the g-width increases 
with the increase of doping. With a much broader mag- 
netic peak in the normal state, the q range (0.35 rlu - 0.65 
rlu) for the normal-state q-scan spectrum (see Fig. 2B of 
Ref. [9]) is too narrow to get meaningful estimates of the 
nonmagnetic background and the normal-state magnetic 
intensity. Moreover, the normal-state g-scan spectrum 
of YBa 2 Cu306.8 is nearly featureless for the same nar- 
row q-range (0.35 rlu - 0.65 rlu), similar to the data of 
Fig. 2B of Ref. [9]. If one would fit the normal-state q- 
scan data only in this narrow g-range (0.35 rlu - 0.65 rlu) 
for YBa 2 Cu306.8 and assume the same g-width as that 
in the superconducting state, one would find that the 
normal-state magnetic intensity is underestimated by a 
factor of 4. Such a significant underestimate may be also 
true for the normal-state magnetic intensities of Tl-2201. 
The authors of Ref. [9] should have extended their mea- 
surements to a wider q range to get reliable estimates of 
the nonmagnetic background and the normal-state mag- 
netic intensities. 

Since magnetic signals in both normal and supercon- 
ducting states strongly depend on doping, as clearly seen 
in YBCO [4], the comparison of the resonance peak inten- 
sities and spectral weights among different compounds 
should be carefully made. For this Tl-2201, there is 
some indication of a broad peak centered around Qaf 
even above T c (Ref. [9]), as observed in underdoped 
YBCO [4,6]. This suggests that the Tl-2201 is slightly 
underdoped so that the magnetic intensity in the nor- 
mal state should be close to that for slightly underdoped 
YBa 2 Cu 3 06.g 2 . This is because both compounds have 
a similar ratio T c /T cm and thus a similar doping level 
(where T cm is the superconducting transition at optimal 
doping). As seen from Fig. la, the normal-state mag- 
netic intensity at 40 meV is about 100 ^/eV per Cu in 
slightly underdoped YBa 2 Cu 3 06.g 2 . The resonance peak 
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intensity for the Tl-2201 can be estimated to be about 
120 /i^/eV per Cu from the resonance spectral weight 
of 0.7fi% per Cu (Ref. [9]). Therefore, the resonance 
peak intensity for the Tl-2201 is close to the normal- 
state magnetic intensity for the slightly underdoped 
YBa2Cu306.92- Since the normal-state magnetic inten- 
sity for Tl-2201 should be similar to that for the slightly 
underdoped YBa 2 Cu 3 6.92, then I s {E r )/I N {E r ) ~ 1.2 
for this single-layer Tl-2201. This is in good agreement 
with feature (e) deduced independently from the differ- 
ence spectrum (Fig. 2a). Moreover, we find that the res- 
onance peak intensity for the Tl-2201 (^120 ^/eV per 
Cu) is a factor of about three smaller than the resonance 
peak intensity for the slightly underdoped YBa2Cu3 06.92 
(^330 H%/eV per Cu). Similarly, the resonance spectral 
weight for the Tl-2201 is also a factor of three smaller 
than that for the slightly underdoped YBa2Cu306.g 2 . 

These important features we have identified above 
should place strong constraints on theories for high- 
temperature superconductivity in cuprates. Any correct 
theories should be able to explain all the magnetic reso- 
nance features in a consistent and quantitative way. In 
a more exotic approach [10], the neutron data are inter- 
preted in terms of a collective mode in the spin-triplet 
particle-particle channel, which couples to the particle- 
hole channel in the superconducting state with d-wave 
OP. This model predicts that the resonance peak energy 
E r is proportional to the doping level p, in disagreement 
with feature (d): E r does not increase with increasing 
p in the overdoped range but is proportional to T c as 
E r /k B T c ~ 5.2 (Ref. [8]). Moreover, this model pre- 
dicts [20] that E r > 2Am (where Am is the maximum 
d-wave gap), which contradicts experiment. Other the- 
ories based on spin-fermion interactions also show that 
E r increases monotonically with increasing p [11,12], in 
disagreement with feature (d). 

Alternatively, feature (d) is consistent with a simple 
particle-hole excitation across the superconducting gap 
within an itinerant magnetism model. This is because 
the particle-hole excitation energy increases with the su- 
perconducting gap which in turn should be proportional 
to T c at least in the overdoped range. This itinerant mag- 
netism model is also supported by very recent Fourier 
transform scanning tunnelling spectroscopic (FT-STS) 
studies on a nearly optimally doped BSCCO [21], which 
show that the quasiparticles in the superconducting state 
exhibit particle-hole mixing similar to that of conven- 
tional Fermi-liquid superconductors. These FT-STS re- 
sults thus provide evidence for a Fermi-liquid behavior in 
the superconducting state of optimally doped and over- 
doped cuprates. Here we quantitatively explain all these 
neutron data [1-9] in terms of an order parameter (OP) 
that has an extended s-wave symmetry [17,18] and oppo- 
site sign in the bonding and antibonding electron bands 
formed within the CU2O4 bilayers [14]. In our model, the 
neutron resonance peak is due to excitations of electrons 
from the extended saddle points below the Fermi level to 



the superconducting gap edge above the Fermi level. 

Within the itinerant magnetism model, neutron scat- 
tering intensity at a wavevector q and an energy E is 
proportional to the imaginary part of the dynamic elec- 
tron spin susceptibility, x" (<fi E) . Qualitatively, with the 
neglect of the Bardccn-Coopcr-Schricffer (BCS) coher- 
ence factor, the bare imaginary part of spin susceptibil- 
ity; Xo (<f> E), i s proportional to the joint density of states 
m E) = E £ S(E E U? E % ), where E % = ^ + A| 
is the quasiparticle dispersion law below T c , is the 
electronic band dispersion, and Ag is the order parame- 
ter [14]. The two particle energy E2(k,q) = E^ + ~ + E^ 
has a minimum and several saddle points for a fixed q. 
The minimum defines the threshold energy, or spin gap 
energy E g , which is achieved at vector k and k + q such 
that both k and k + q belong to the Fermi surface and 
thus E g = 2A £ (Ref. [14]). 

The joint density of states has divergences at the ex- 
tended saddle points [14]. The divergent peak in x" (g, w ) 
occurs because of transitions between the occupied states 
located in the extended saddle points below Ep and 
empty quasiparticle states at the superconducting gap 
edge above Ep. The saddle points produce Van Hove sin- 
gularities in the quasiparticle density of states in the su- 
perconducting state at an energy of \J( e ]* H ) 2 + A| and 
the superconducting condensate creates a sharp coher- 
ence peak in the quasiparticle density of states at the 
gap edge. Thus, the divergence in the joint density of 
states in the superconducting state is then located at the 
energy E* = A^ + _.+ [eX H ) 2 + A|. This simple expres- 
sion for E* has been verified by numerical calculations in 
the case of an isotropic s-wavc order parameter [14]. 

In Fig. 3, we plot the Fermi surface for a slightly under- 
doped BSCCO, which is inferred by extending a portion 
of the Fermi surface determined by angle-resolved pho- 
toemission spectroscopy (ARPES) [22]. There are ex- 
tended saddle points that are located near (±7r, 0) and 
(0, ±tt), as shown by ARPES [23]. A solid line in Fig. 3 
represents a segment of the extended saddle points very 
close to the Fermi surface. Other portions of the saddle 
points [e.g., along the line from (0, 0) to (n, 0) ] are signif- 
icantly away from the Fermi surface [23] and not shown in 
the figure. If we only consider the magnetic excitations at 
a fixed wavevector that corresponds to the antiferromag- 
netic wavevector Qaf, only four electron wavevectors at 
the Fermi surface are connected by Qaf as indicated by 
arrow 1 in Fig. 3. Each of these vectors forms an angle of 
r with respect to the Cu-0 bonding direction. The elec- 
tron transitions from the occupied states located at the 
extended saddle points below Ep to empty quasiparticle 
states at the superconducting gap edge above Ep are in- 
dicated by arrow 2. Because the quasiparticle densities 
of states at the gap edge and the saddle points are diver- 
gent at zero temperature, such transitions will produce a 
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sharp resonance peak at E r 
(Ref. [14]). 
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FIG. 3. The Fermi surface for a slightly underdoped 
BSCCO with T c = 88 K. This Fermi-surface is extrapolated 
from a part of the Fermi surface determined by ARPES [22] 
using symmetry arguments. Arrow 1 indicates electron tran- 
sitions from the occupied states in the superconducting gap 
edge below Ef to empty quasiparticle states at the gap edge 
above Ef. Arrow 2 marks electron transitions from the occu- 
pied states located in the extended saddle points below Ef to 
empty quasiparticle states in the superconducting gap edge 
above Ef- 

In terms of 9 r , both E 9 (Qaf) and E t {Qaf) can be 
rewritten as 



E g {Q AF ) - 2A(9 r ) 



and 



E r (Q AF ) = A(9 r ) + yJ[A(0 r )]* + (t 



Here y/[A(9 r )] 2 + (eY H ) 2 is the energy of a saddle point 
below the Fermi level along the 9 r direction. One should 
note that Eq. 2 is valid only if the saddle points are very 
close to the Fermi surface, as is the case. 

We now consider the BCS coherence factor that has 
been ignored in the above discussions. The BCS coher- 
ence factor is given by [24] : 



e k £ k+q~ 



^k+q^k \ 



E k E k+q 



For ££<< A^ and zj; +q - « tne BCS coherence 

factor is close to 1 when Ar + ~ and A? have opposite 
sign, but is close to zero when Ar + ~ and A^ have the 
same sign. The coherence factor is not negligible for 
~ Ar and er , - = even if A.- . - and Af have the same 

k k+q k+q k 

sign. 

For a single-layer compound with d-wave order pa- 
rameter symmetry, Ar + ^ Ap and Ar have opposite 



./(eY H ) 2 + Al sign so that the BCS coherence factor ~ 1 and thus 
V k k I s (E r )/I N (E r )» 1 (see Ref. [13]). Therefore, the ex- 

perimental observation of I s (E r ) j I N (E r ) ~ 1.0 in the 
single-layer Tl-2201 [9] rules out the d-wave order pa- 
rameter. Alternatively, for a single-layer compound with 
an s-wave symmetry, A-r + ^ Ap and A^ have the same 
sign, so the BCS coherence factor could be much less 
than 1. This may lead to I s (E r ) / I N (E r )~ 1, in agree- 
ment with feature (e). Hence, only the intralayer (intra- 
band) s-wave symmetry is compatible with feature (e): 
I s (E r ) j I N (E r ) ~ 1.0. On the other hand, if extended 
saddle points are significantly below the Fermi level, the 
BCS coherence factor will be substantial even for the in- 
tralayer (intraband) magnetic scattering in the case of 
s-wave gap symmetry. Then I s (E r )/I N (E r ) could be 
much larger than 1 in this special case. Therefore, the 
observation of I s (E r )/I N (E r ) >> 1 for the intralayer 
(intraband) magnetic scattering is consistent with either 
s-wave or d-wave gap symmetry. 

For a double-layer compound, interactions within 
CU2O4 bilayers yield bonding and antibonding bands. 
Transitions between electronic states of the same type 
(bonding-to-bonding or antibonding-to-antibonding) and 
those of opposite types are characterized by even or odd 
symmetry, respectively, under exchange of two adjacent 
Cu02 layers. As a result, the magnetic excitations be- 
tween different bands correspond to odd channel excita- 
tions while the excitations within the same band to even 
channel excitations [4]. If the order parameters in the 
bonding and antibonding bands have the same sign, the 
magnetic resonance intensities in both channels should 
be similar for an intraband pairing symmetry of cither 
s-wave or d-wave. This is because the BCS coherence 
factors for both channels are the same in this case. On 
the other hand, if the order parameter has an s-wave sym- 
metry and opposite sign in the bonding and antibonding 
electron bands, the magnetic resonance intensity in the 
odd channel will be much larger than that in the even 
channel due to the large difference in their BCS coher- 
ence factors. Therefore, only if the order parameter has 
an s-wave symmetry and opposite sign in the bonding 
and antibonding electron bands, can the observed feature 
( c ) : I odd( E r)/ I even( E r) = 5 for the slightly underdoped 
YBCO and I^ dd (E r )/I^ ven (E r ) >10 for the overdoped 
YBCO [3] be explained within the itinerant magnetism 
approach. 

Based on the t-J model and d-wave pairing symmetry, 
Brinckmann and Lee [25] have recently attempted to ac- 
count for feature (c) using a very unrealistic parameter: 
J±/J\\ = 0.6 (where J\\ and Jj_ are the effective intralayer 
and interlayer antiferromagnetic exchange energies, re- 
spectively). For undoped YBCO, neutron experiments 
[4] show that J±/J\\ = 0.1, which is far less than 0.6 
used in the calculation [25]. Further, the value of J±/J\\ 
will be negligible if Jn remains substantial and the opti- 
cal magnon gap (oc y/J\\J±) goes to zero, which should 
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be the case for optimally doped and overdoped YBCO. 
Neutron data for underdopcd YBCO [26] show that the 
optical magnon gap is about 50 meV for YBa 2 Cu 3 6 .5 
and is reduced to about 25 meV for YBa 2 Cu 3 06.7- If 
we linearly extrapolate the optical magnon gap with the 
oxygen content, the gap will tend to zero in YBa 2 Cu30 y 
for y > 0.9, implying that J±/J\\ « 0.1 for optimally 
doped and overdoped cuprates. 

Using a more realistic parameter of J± and the d-wave 
pairing symmetry, Millis and Monien [27] appear to be 
able to explain feature (c). They showed that the reso- 
nance spectral weights in the odd and even channels are 
related to the difference between the resonance energy 
E r and the particle-hole spin excitation energy 2A(6> r ) 
at Q AFl i.e., w odd /W even = [2A(6 r ) - E° dd ]/[2A{6 r ) - 
E™ en ]. Here 2A(0 r ) must be larger than both E% ven and 
E odd (Rcf. [27]). The INS data of Yo.gCao.iBaaCugOr-j, 
would be consistent with this theoretical model if one 
would choose an unrealistic parameter 2A(8 r ) = 49 meV 
(Rcf. [19]). As discussed below, 9 r is about 15° for 
slightly overdoped cuprates so that 2A(8 r ) =1.73Am 
(where Am is the maximum d-wavc gap). The intrinsic 
tunneling spectra, which are unsusceptible to surface de- 
terioration, can provide the most reliable determination 
of the bulk superconducting gap [28] . From the intrinsic 
tunneling spectra, one finds that 2Am /k B T c = 8.09 for 
the optimally doped BSCCO with T c = 94 K (Ref. [28]), 
2A M /k B T c = 6.73 for a slightly overdoped BSCCO with 



T c = 89 K (Ref. [28]), and 2A M /k B T c = 5.37 for an over- 
doped BSCCO with T c = 80 K (Rcf. [29]). It is apparent 
that 2A M /k B T c decreases almost linearly with T c in the 
overdoped range. Using the fitted curve of 2Aulk B T c 
versus T c , we estimate 2Am /k B T c = 6.04 and 2A(9 r ) = 
38.6 meV for Y .9Ca .iBa 2 Cu 3 O7-2 / with T c = 85.5 K. 
Similarly we can obtain 2Am /k B T c = 5.66 and 2A(0 r ) = 
35.1 meV for overdoped BSCCO with T c = 83 K. INS ex- 
periment on Yo.9Cao.iBa 2 Cu307_j / (Ref. [19]) indicates 
E even = 43 mcV > 2A(6 r ), in contradiction with the 
theoretical model [27]. For overdoped BSCCO with T c 
= 83 K, E odd = 38 mcV (Ref. [8]) and thus E^ ven = 45 
meV by analogy with the case of Yo.9Cao.iBa 2 Cu307_ !r 
Both Ef en and E° dd in this overdoped BSCCO arc far 
larger than 2A(# r ), in contradiction with any theoretical 
models based on the d-wave pairing symmetry. 

For the underdopcd YBa2Cu306.7, the even-channel 
magnetic intensity in the normal state is a factor of 2.0- 
2.5 lower than the odd-channel one for E = 40 meV, 
which is about 15 meV above the optical magnon gap 
[26]. This implies that the intcrlayer antiferromagnetic 
correlation does not influence magnetic excitations well 
above the optical magnon gap [26]. Since the opti- 
cal magnon gaps for optimally doped and overdoped 
YBCO arc much smaller than that for the underdopcd 
YBa 2 Cu 3 06.7, one should expect that I? M (E) / ' I* en {E) 
~ 1 for E ~ 40 meV in optimally doped and overdoped 



TABLE I. Comparison of experiments with extended s-wave, d-wave, and isotopic s-wave. In both extended and isotropic 
s-wave models, the order parameters are assumed to have opposite signs in the bonding and antibonding electron bands formed 
within the CU2O4 bilayers. Here DA = definitive agreement, A = agreement, QA = qualitative agreement, PA = possible 
agreement, D = disagreement and DD = definitive disagreement. 



INS data 


extended s-wave 


d-wave 


isotropic s-wave 


Feature (a): I s (E r ) / 1" (E r ) = 3.6 for YBCO 


DA 


DA 


DA 


Feature (c): I^ dd (E r )/I^ en (E r ) = 5 for YBCO 


DA 


DD 


DA 


Feature (e): I s (E r )/ I N (E r ) <1 for Tl-2201 


DA 


DD 


DA 


The magnitudes of E r and E g 


DA 


DD 


DD 


Spin gap in La2-xSr x Cu04 


DA 


DD 


DD 


Other data 


extended s-wave 


d-wave 


isotropic s-wave 


Penetration depth 


DA 


QA 


DD 


Thermal conductivity 


DA 


QA 


DD 


Specific heat 


DA 


QA 


DD 


Nonlinear Meissner effect 


DA 


QA 


DD 


ARPES (underdoped) 
ARPES (overdoped) 


A 

DA 


A 

DD 


DD 
DD 


Quasiparticle tunneling (underdoped) 
Quasiparticle tunneling (overdoped) 
Raman scattering (underdoped) 


PA 
DA 
PA 


PA 
DD 
PA 


DD 
DD 
DD 


Raman scattering (overdoped) 


DA 


DD 


PA 


NMR/NQR 
Andreev reflection 


A 
A 


A 
PA 


DD 
DD 


Pb/c-axis YBCO Josephson junction 


A 


D 


A 


c-axis BSCCO twist Josephson junction 


DA 


DD 


PA 


Corner SQUID/ Josephson junction 
Tricrystal/Tetracrystal Josephson junctions 


PA 
PA 


PA 
PA 


PA 
PA 
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YBCO. Thus the observed feature 
(c): Iodd( E r)/Ieven(Er) >> 1 for optimally doped and 
overdoped YBCO can only be explained by an OP that 
has s-wave symmetry and opposite sign in the bonding 
and antibonding electron bands. 

For the slightly undcrdoped YBa 2 Cu3 06.92, we can 
also deduce the value of I^ ven (E r ) / I^ ven (E r ) using the 
measured I^ dd (E r )/I^ dd (E r ) = 3.6, I^ dd (E r )/I^ ven (E r ) 
= 5, and the inferred lodd(-^r) / I^ ven {E r ) — 1 (see above 
discussion). From the measured Ig dd {E r ) / I^ ven {E r ) 
= 5, we have I^ ven (E r ) = Q.2I° dd {E r ) and thus 
lLn(E r )/lLn(E r ) * Q.2I s odd {E r )/I? dd {E r ) = 0.72. 
This value is close to that deduced for the single-layer 
Tl-2201 (~ 1). These results consistently suggest that 
the intraband neutron intensity at E r in the supercon- 
ducting state is close to that in the normal state. This 
unique and important feature we have identified rules out 
d-wave order parameter symmetry because d-wave sym- 
metry predicts [13,25] that lf ven {E r ) / I^ ven (E r )»l for 
bilayer compounds and I s (E r )/I N (E r )»l for single- 
layer compounds. 

From Eqs. 1 and 2, it is easy to calculate E g and E r if 
one knows the gap function A(8) and the 8 r value. From 
the measured Fermi surface, one can readily determine 
8 r . For example, we find 8 r = 18.4° for a slightly un- 
dcrdoped BSCCO from Fig. 3. For optimally doped 
cuprates, we get 8 r ~ 16.0°. If we would use an isotropic 
s-wave gap function A(8) = 28 mcV for a slightly over- 
doped YBCO, we would have E g = 56 meV and E r > 
56 meV, which are far larger than the measured E g = 
33 mcV and E r = 40 meV (Ref. [3]). If we would use 
a d-wave gap function A (8)— 28cos2(9 meV, we would 
have Eg = 47.5 meV and E r > 47.5 mcV, which are 
also far larger than the measured values. Thus, one can- 
not quantitatively explain the neutron data in terms of 
d-wave and isotropic s-wave symmetries. 

Alternatively, an extended s-wave with eight line nodes 
(Ai g symmetry) is in quantitative agreement with two- 
thirds of the experiments that were designed to test the 
order-parameter symmetry for hole-doped cuprates [17]. 
The remaining one-third (e.g., tricrystal grain-boundary 
Josephson junction experiments) are explained qualita- 
tively by Zhao [17] and by Brandow [18]. In Table 1, 
we compare nearly all the experiments used to test the 
OP symmetry with the extended s-wave, d-wave, and 
isotropic s-wave models. In both extended s-wave and 
isotropic s-wave models, the order parameters are as- 
sumed to have opposite sign in the bonding and anti- 
bonding electron bands formed within the CU2O4 bilayers 
[14]. The detailed comparisons with other experiments 
are made in Ref. [17]. From Table 1, one can see that the 
neutron data alone provide a definitive answer to the in- 
trinsic, bulk OP symmetry because INS is a bulk, phase 
and angle sensitive technique. Other bulk and non-phase 
sensitive experiments provide complementary support to 
the present conclusions. The phase and surface sensitive 
experiments cannot definitively determine the intrinsic 



bulk OP symmetry because the surface OP might be dif- 
ferent from the bulk one [17]. 

For a slightly overdoped YBCO, more than six inde- 
pendent experiments consistently suggest that [17] the 
gap function is A(8) = 24.5(cos46> + 0.225) mcV. Sub- 
stituting 8 r — 16° into the gap function, we get A(8 r ) 
= 16.3 mcV and thus E g = 32.6 meV, in quantitative 
agreement with the measured one (32-33 meV). 

In our model, the position of the extended saddle point 
along the 8 r direction is located at E r — E g /2 in the su- 
perconducting state (see Eqs. 1 and 2). For optimally 
doped YBCO with E r = 41 meV and E g = 32 mcV, we 
find that the saddle point in the superconducting state 
is located at an energy of 25 meV below the Fermi level. 
This is in agreement with the ARPES studies [30] which 
suggest that the Fermi level in the superconducting state 
for optimally doped cuprates is < 30 meV above the ex- 
tended saddle points that have the same energy over a 
large momentum space [23]. Further, electronic Raman 
scattering spectra in YBa2Cu40s have been used to de- 
termine the energy of the extended saddle points more 
accurately [31]. At 10 K (well below T c ), the energy of 
the extended saddle points is found to be 24.3 meV be- 
low the Fermi level, in excellent agreement with the result 
predicted by our model from the INS data (~ 25 meV). 

We have identified the unique and important feature: 
I (E r )/I N (E r ) ~ 1.0 for the intralayer and intraband 
magnetic scattering. This feature unambiguously rules 
out the d-wave OP symmetry. The unambiguous deter- 
mination of the intrinsic extended s-wave pairing symme- 
try for hole-doped cuprates places strong constraints on 
the pairing mechanism for high-temperature supercon- 
ductivity. A recent calculation suggests that high energy 
Cu-0 charge fluctuations can lead to an attractive in- 
teraction between conduction electrons and the pairing 
symmetry may be of extended s-wave (A\ g ) [32]. 

Indeed, precise thermal-difference reflectance spectra 
of several cuprate superconductors (T c =105-120 K) ex- 
hibit pronounced features at photon energies of about 2.0 
cV, which may be related to the Cu-0 charge fluctuations 
(Ref. [33]). These features can be well described within 
Eliashberg theory with an electron-boson coupling con- 
stant X c h of about 0.40. In order to explain a supercon- 
ducting transition temperature of 105-120 K, the authors 
simulated [33] an electron-phonon coupling feature at 50 
meV with the coupling constant X p h of about 1.0. This 
Eliashberg model's simulated electron-phonon coupling 
agrees well with the results obtained from first principle 
calculations [34]. 

The contribution of this 2 eV component can be equiv- 
alent to a negative Coulomb pseudopotential \i* within 
Eliashberg theory [35]. By using a realistic electron- 
phonon coupling spectral weight deduced from tunneling 
spectra and a [i* = -0.15, and taking into a polaronic 
effect, Zhao et al. [35] are able to explain the negligible 
isotope effect on T c , and the magnitudes of T c and the 
superconducting gap for optimally doped 90 K supercon- 
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ductors. 

In summary, we have identified several important fea- 
tures in the neutron scattering data of cuprates, which 
arc difficult to be explained in terms of d-wave and 
isotropic s-wave order parameters. Alternatively, we 
show that the neutron data are in quantitative agreement 
with an order parameter that has an extended s-wave 
(Ai g ) symmetry and opposite sign in the bonding and an- 
tibonding electron bands formed within the CU2O4 bilay- 
ers. This A lg pairing symmetry may be compatible with 
a charge fluctuation mediated pairing mechanism. High- 
temperature superconductivity in cuprates may be due to 
the combination of strong electron-phonon coupling and 
substantial electron-charge-fluctuation coupling [33,35]. 

* Correspondence should be addressed to 
gzhao2@calstatela.edu 
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